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\ Within the harmonic approximation, the effects of lattice vibration on the thermodynamics of 
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^ I nano-sized coherent clusters in bcc-Fe consisting of vacancies and/or copper are investigated. A 

■ combination of on-lattice simulated annealing based on Metropolis Monte Carlo simulations and off- 

I . lattice relaxation by Molecular Dynamics is applied to obtain the most stable cluster configurations 

^1 at T = K. The most recent interatomic potential built within the framework of the embedded 

1^1 ' atom method for the Fe-Cu system is used. The vibrational part of the total free energy of defect 

. clusters in bcc-Fe is calculated using their phonon density of states. The total free energy of pure 
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, bcc-Fe and fcc-Cu as well as the total formation free energy and the total binding free energy of 

I the vacancy-copper clusters are determined for finite temperatures. Our results are compared with 

• the available data from previous investigations performed using empirical many-body interatomic 

^—1 . potentials and first-principle methods. For further applications in rate theory and object kinetic 
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^ I Monte Carlo simulations, the vibrational effects evaluated in the present study are included in the 

^ ■ previously derived analytical fits based on the classical capillary model. 
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I. INTRODUCTION 



The level of impurity copper in reactor pressure vessel (RPV) steels of existing fission 



reactors (0.03 to 0.3 wt% 

n 



n 



is typically higher than the solid solubility of Cu in Fe at RPV 



operation temperature [l|-|3|. Neutron irradiation gives rise to a steady-state vacancy con- 
centration that is orders of magnitude higher than in thermal equilibrium and enables Cu 
to diffuse efficiently via a vacancy mechanism 0, 5|. As a consequence, Cu-rich clusters 
form and take up a number of vacancies. In fact, mixed Cu- vacancy clusters with non- 
monotonically varying fractions of Cu and vacancies were observed recently in binary Fe-Cu 
alloys js Sl. Similar processes take place in RPV steels, although other impurity and alloy- 
ing elements may also play a role and vacancies seem to contribute less. Cu-rich clusters 
or precipitates in both binary Fe-Cu alloys (Cu <0.3 wt%) and RPV steels are typically 



smaller than 4 nm in diameter 



10| . Experimental investigations UJ, ll2[ also show that 



these Cu-rich precipitates are coherent i.e. they possess the bcc-structure of iron. 

Irradiation-enhanced formation and evolution of the Cu-rich nano clusters or precipitates 
are multiscale processes. These processes involve various physical phenomena and can be 
efficiently investigated with rate theory [s], [isl and object kinetic Monte Carlo (OKMC) 
simulations [l^ . Iisl covering more or less realistic time and length scales. Important ther- 
modynamic parameters, in particular the binding free energy of monomers (Cu-atoms or 
vacancies) to defect clusters required for these coarse-grained or continuum approaches are, 
however, hardly accessible from experiments. But atomistic computer simulations can pro- 
vide these parameters. 

In this study, our main objective is to quantitatively determine the vibrational contribu- 
tion to the thermodynamics of defect clusters consisting of vacancies and/or Cu-atoms in 
the Fe-Cu system with cluster size larger than that first principle methods can treat. For 
this purpose, the most recent interatomic potential for the Fe-Cu alloy is employed 16|. We 
perform here a combination of Metropolis Monte Carlo (MMC) simulations and Molecu- 
lar Dynamics (MD) simulations to find the energetically most stable defect configurations. 
This method is previously employed by Al-Motasem et al. 17| and is similar to the proce- 
dures used by Takahashi et al. 18[ and Kulikov et al. 19|]. For small Cu-clusters (up to 
15 Cu-atoms) in bcc-Fe, first principle studies using density functional theory (DFT) reveal 
that the vibrational effect is comparable (~36%) to configurational entropy |20|. To deter- 



mine the vibrational effects, we find the phonon frequencies of the stable structures under 
the harmonic approximation, with the computation of thermodynamic quantities following. 
The total formation free energy and the total binding free energy for pure and mixed defect 
clusters are calculated. Moreover, in order to furnish the input parameters for further ap- 
plications in rate theory and OKMC simulation, we also aim here to determine simple but 
practically correct analytical fit formulae that render the total and the monomer binding 
free energy of defect clusters for any finite temperatures. The computational procedure 
is detailed in section 2. Results are presented and discussed in section 3, followed by the 
summary and conclusion of this work. 



II. COMPUTATIONAL METHOD 



In order to find the most stable configurations of coherent clusters consisting of I vacancies 
and m Cu atoms (i.e. viCum clusters), at first we performed on-lattice simulated annealing 
based on Metropolis Monte Carlo (MMC) simulation as described in ref. The initial 
configurations contained randomly-distributed isolated vacancies and Cu atoms at 600 K. 
Then, we decreased the temperature to K reducing it by 60 K at every 100 MMC steps. 
A single cluster is found at K for all cases considered in this study. Subsequently, we 
performed a quasi-dynamic quench through Molecular Dynamics simulation at K to relax 
the positions of Fe and Cu atoms off-lattice following the procedure reported in ref. [l^ . The 
lowest formation energy, Ef{viCum) is then obtained by performing several such relaxations 
for shghtly different values of the lattice constant. This quantity is defined as 

Ef{viCum) = E{viCum) - {N - I - m)E'j,f - mE^"^ (1) 

where is the total number of regular lattice sites. E{viCum) is the configurational energy 
of the simulation box containing the most stable configurations of the viCum cluster. E'^^ 
and i?™^ represent the cohesive energy per atom of bcc-Fe and fcc-Cu, respectively {E'^p^ 
= -4.122 eV and E^^ = -3.540 eV for chosen interatomic potential) which are chosen as 
reference systems. The total binding energy of a viCum cluster, Eund{viCum) describes the 
energy change when isolated vacancies and Cu atoms are combined to form a cluster and is 
defined by 

Eund{VlCUm) = Ef{viCUm) " lEf{vi) - mEf{Cui), (2) 



where Ef{vi) and Ef{Cui) are the formation energy of a monovacancy and of a single Cu 
atoms, respectively (cf. Eq. ([1])). Ef{Cui) is also called substitutional energy. Negative 
values of Ei^indiyiCum) suggest an energetically favorable condition for the monomers to form 
a cluster. We employed here, the most recent and currently the rnost suitable empirical 
many-body interatomic potential of Pasianot and Malerba (PM) [161] developed for the 



description of radiation effects in Fe-Cu binary alloys. It is based on the Mendelev 2ll] and 



the Mishin [2j] potentials for pure Fe and Cu, respectively. The simulation box of bcc-Fe 
lattice is chosen to be a cubic in shape with each edge of 10a length, where a is the lattice 
constant. Three dimensional periodic boundary conditions are applied. In this study, pure 
clusters containing up to 80 either vacancies or copper atoms are considered. Maximum 
cluster size amounts to total of 100 vacancies and copper atoms for mixed clusters. In 



agreement with experimental findings 



11 



12| . all the clusters considered in this paper are 



assumed to possess the bcc-structure of iron. 

The phonon frequencies Ui are determined by diagonalizing the dynamical matrix under 
the harmonic approximation for the most stable configurations obtained at K. At F-point, 
the dynamical matrix is defined as 

*^. = -(^^f^)> (3) 

where rrii is the mass of z*^ atom. It is calculated by applying small displacements rj away 
from their most stable atomic positions and evaluating the subsequent induced forces (Fi) 
and force derivatives (^). 

The phonon calculations allow us to compute various thermodynamic properties, e.g. 
free energy, entropy, heat capacity etc. The vibrational contribution to the total free energy 
(from now on the vibrational free energy) of a solid, at finite temperature T is calculated 
using the following expression: 

Qvib^y^ T) = U"\V, T) - TS^'^V, T), (4) 

where ?7™^(V,T) represents the vibrational internal energy and S^'^V^T) the vibrational 
entropy. In the harmonic approximation, these contributions (cf. Eqs. ElE])) are represented 
as the sum of individual normal frequencies, ojj. The total number of individual normal 
frequencies are three times the total number of atoms in the primitive cells. Because of the 
translational invariance of the system three eigen values corresponding to the translational 
degrees of freedom are zero and hence, 



and 



3Af-3 



gfegT — 1 _ 1 — e ''bt 



(6) 



4 = 1 

Here, ks represents the Boltzmann constant. 

An alternative way to calculate the vibrational free energy [G""^^) here, is through the 
phonon density of states (PDOS), g{u): 

T) = (3iV - 3) ksT In (^2stnh ) 9ioj)duj, (7) 

where g{oj) is the normalized so that 

POO 

/ giu)du = 1. (8) 

The phonon density of states g{oj), which can also be considered as a histogram of phonon 
frequencies with vanishing bin size, gives the number of vibrational modes per phonon 
frequency. 

The vibrational part of total formation free energy of a viCum cluster is calculated from 
the following expression: 

Gf{viCum. V, T) = G-'\viCum. v,T)-{N-l- m)GZi{Fe, V, T) - niGZHCu, V, T), (9) 

where, G^^^{viCum,V,T) determined using Eq. ([7]) represents the vibrational free energy 
of the system with the viCum cluster. (^^^^(Fe, T) and G^^J'f^{Cu,V,T) represent the 
vibrational contribution to the cohesive energy per atom at given temperature (T) and 
volume (V) of bcc-Fe and fcc-Cu, respectively. The total formation free enery is then written 

as 

G'™(t;;CM^, V, T) = EfiviGum) - TSj'^f {viGum) + Gf{viGum, V, T), (10) 

where Sj'"'^ (= kslnW) is the formation entropy at T = K which is determined by the 
number W of different three-dimensional arrangements of a cluster with a given shape at 
a given position. In the present paper the contribution of Sj""'^ {viCum) is neglected since 



the computation time is prohibitively long. The efficient determination of Sj'"' {viCum) is 
a subject of ongoing investigation. 

The vibrational part of total binding free energy and the total binding free energy can 
be obtained in a similar fashion from the following formulae, 

G'ilUviCurn, V, T) = Gf{viCu^, V, T) - IGj'iv,, V, T) - mGf{Cu,, V, T) (11) 

and 

Gltnd\viCu^, V, T) = EMndiviGum) + Gltdi^iGu^, V, T) (12) 

where Gf\viGu^,V,T) , Gf\vi,V,T) and Gf''{Gui,V,T) represent the vibrational part 
of free formation energy of the ViCum cluster, of a monovacancy and a single Cu in bcc-Fe, 
respectively (cf. Eq. f lTO|) ). Eund{viGum) represents the binding energy of a cluster at K 
(cf. Eq. dlD). 

The monomer binding free energy is a binding energy of a monovacancy or of a single 
Cu atom to a pure vi or Gu^ cluster. This quantity, at T 7^ and zero pressure is required 
as input parameters for rate theory and object kinetic Monte Carlo simulations and can be 
determined from the following expression: 



Gl°"^\Gu, Cu^, V, T) = G™(C^z^, V, T) - Gl^^fiGu^.,, V, T) (13) 



III. RESULTS AND DISCUSSION 



A. Phonon density of states and vibrational free energy of bcc-Fe and fcc-Cu 



Figure 1 depicts the phonon density of states (PDOS) of pure bcc-Fe and fcc-Cu, which 
are in good agreement with the available experimental observation and theoretical calcula- 



tions 



23 



26j. The PM potentia 



and on the Mishin potential 



which is based on the Mendelev potential j2lj for pure Fe 
22| for pure Cu seems to reproduce the experimental phonon 
density of states of elemental bcc-Fe and fcc-Cu well. The vibrational free energy of the 
considered systems is then calculated using their PDOS in Eq. ([7]). As shown in Fig. 2, our 
results for bcc-Fe and fcc-Cu are compared with the Scientific Group Thermodata Europe 
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FIG. 1. 

Phonon density of states of (a) bcc-Fe and (b) fcc-Cu. 



(SGTE) database and results from other empirical potentials for these two elements 27- 
3l[ |. The SGTE methodology includes CALPHAD calculations and experimental data and 



provides the thermodynamic data for inorganic and metallurgical systems. The CALPHAD 
data in re f. |28| has been obtained with the THERMOCALC program and the SGTE unary 
database 32]. The vibrational free energy of bcc-Fe and fcc-Cu calculated in the present 
work fall very much into the range of available literature data. At 1200 K, the calculated vi- 
brational free energy of bcc-Fe using the PM potential is about 160 meV and 100 meV higher 
than SGTE database [^^j and CALPHAD data 28], however, it falls within 50 meV range of 



vibrational free energy calculated using other empirical potentials such as the Meyer-Entel 
potential, the Johnson potential which is a recent update of the Johnson-Oh potential, and 
the Finnis-Sinclair potential At around 600K - an operating temperature of a nuclear 
reactor, the former difference in the vibrational free energy is about 70 meV and 30 meV, 
respectively while the latter is within 20 meV range. The Meyer-Entel potential predicts the 
highest values of the vibrational free energy and hence the farthest from SGTE database or 
CALPHAD data while the PM potential provides the closest values. Similarly, though the 



Ludwig-Farkas potential predicts the vibrational free energy of the fcc-Cu (in 30|]) closest 



to the SGTE database 23], the PM potential calculates it within the closer range of 30 meV 
and 10 meV at 1200 K and 600 K respectively. The ability to reproduce the experimen- 
tal PDOS of bcc-Fe and fcc-Cu to a satisfactorily extent and to predict their vibrational 
free energy close to the SGTE database, the choice of the PM potential is justified for the 



considered systems in the present work. 




FIG. 2. Vibrational free energy (meV/atom) of (a) bcc-Fe and (b) fcc-Cu. Present results with 
the PM potential are compared with the SGTE or CALPHAD data and the other potentials for 



bcc-Fe 



27 



29( 1 and for fcc-Cu 



B. Formation free energy of the clusters 

The presence of considered number of vacancies and substitutional Cu-atoms seems to 
very weakly affect peak-positions and form of the PDOS of bcc-Fe to notice any visible 
difference. For this reason, the PDOS of defect clusters in bcc-Fe are not shown here. First 
principle study for a monovacancy in bcc-Fe also reports a similar observation with a slight 



phonon spectrum-shift towards lower frequencies [33|. Figure 3 depicts the total formation 
free energy of a monovacancy and a single Cu atom in bcc-Fe. In the present study, the 
values of formation energy of a monovacancy and of a single Cu atom are Ef{vi) = 1.710 
eV and Ef{Cui) = 0.4369 eV, respectively which are obtained using the PM potential jl?! . 
Due to an increase in vibrational interactions (or thermal excitations) with temperature, the 



probability to form vacancies in a lattice also increases. In other words, the vibrational con- 

7 



tributions lower the total formation free energy of a monovacancy {G'^"^"'''{vi, T)) and of pure 



vacancy clusters T)) at higher temperatures as shown in Fig. 3(a) and Fig. 5(a), 

respectively. Recent computational investigations on the monovacancy in bcc-Fe with both 



first principle calculations 33| and many-body empiricial potential 3J] also conform to our 



results for a monovacancy case, and hence for pure vacancy clusters {vi) in general. However, 



it is worth to note that the formation energy of a monovacancy in bcc-Fe estimated from 



experimental investigations 



35| and using other empirical potentials [34 



different from the first principle calculations 



33 



39 




are slightly 



44j . Mendelev et al. [341 estimated and 



interpolated the temperature dependence of total formation free energy for a monovacancy 
in bcc-Fe. They fitted interpolation coefficients to MD data at low temperatures and the 
equilibrium defect concentrations at high temperatures. These coefficients seem to capture 
additional anharmonic effects leading to a further decrease of G'J°*"'(t'i, T) above 900 K as 
compared to our results. Nevertheless, in agreement with the available data our results for 
the monovacancy case provided us confidence to further the similar computation for pure 
vacancy clusters (vi) of larger size. The decrease in total formation free energy of pure vi 
clusters, demonstrated in Fig. 5(a) as a function of temperature, indicates an enhanced 
activity of vacancy cluster formation or vacancy trapping in bcc-Fe at higher temperatures. 




FIG. 3. Total formation free energy (eV) of (a) a monovacancy in bcc-Fe and (b) a single Cu atom 
in bcc-Fe. 



Our results for total formation free energy of a single Cu-atom in bcc-Fe, on the other 
hand is contrary to our expectation and the prediction of first principle calculations [26| (see 
Fig. 3(b)). First principle calculation using the VASP code predicts the formation energy 



of a single Cu in bcc-Fe, Ef{Cui) = 0.55-0.77 eV [26|, |40|, l4l|, |45|] while the SIESTA code 
predicts Ef{Cui) = 0.48 eV [43|. The rigid lattice model |43[ estimates Ef{Cui) = 0.545 



eV. The PM potential fitted against the experimental values of Cu solubility in bcc-Fe [sl, I^gJ 
and neglecting the vibrational contribution estimates the lower formation energy, Ef{Cui) 



0.4369 eV as compared to the reported DFT values 



40 



41 



45|. In Fig. 3(b) we 



compare the temperature dependence of the total formation free energy of a single Cu-atom 



26|. The values of vibrational part 



in bcc-Fe with the VASP calculations of Reith et al. 
of total formation free energy {Gy''{Cui,V,T)) as reported by Reith et al. are negative 



and the absolute values of Gj''(Cui,V,T) increase with temperatures. As shown in Fig. 
3(b) the total formation free energy, Gj°*'''(Cwi, y, T) remains positive and decreases with 
temperature where the relatively high formation energy [EfiCui) = 0.77 eV) is summed with 
negative values of vibrational part of total formation free energy (G™*(Cmi, V, T)) according 
to Eq. ( ITO|) . The PM potential, on the other hand renders the positive values of Gj**(Ciii, T) 
with its absolute values increasing with temperature. The calculated total formation free 
energy as demonstrated in Fig. 3(b) increases with temperature unlike the case reported 
by Reith et al. Although the increase in total formation free energy is small (~ 0.05 eV at 
1200 K), it indicates a slight unfavourable condition for the formation a single Cu-atom in 



bcc-Fe at higher temperatures. A full thermodynamic integration j47|, |48l| being consistent 
with the harmonic approximation in the present study, also qualitatively predicts a similar 
trend of total formation free energy of a single Cu in bcc-Fe as a function of temperature. 

The solubility limit of Cu in iron or Cu concentration (in at %) in iron at equilibrium, 
C^°' is related to the total formation free energy by 

G'™(C«i,V^,T) = -kBTln{C^"J). (14) 

Figure 4 shows experimental and theoretical solubility data from the literature and they 
are compared with the results obtained from Eq. ( !T4|) . Experimental data were obtained 
from measuring the concentration and diffusion profile of a thin Cu layer deposited on an 



iron crystal 49|; thermoelectric power and small angle X-ray scattering measurements in 
thermally aged pure Fe-Cu alloys jsj. The computational data were obtained using an ad- 
vanced Monte Carlo method implemented in the Alloy Theoretic Automated Toolkit package 
(ATAT) 50| with the PM potential denoted as COS. 20 16|. Note that in both experimental 
and theoretical investigations the solubility was determined for the thermodynamic phases, 
fcc-Cu and bcc-Fe in equilibrium. Our results for solubility limits of Cu in bcc-Fe with 
the PM potential neglecting the vibrational effects i.e. corresponding to the calculations at 
T = K formation energy are in a good agreement with experimental and computational 
[igI, 49 1 . Inclusion of vibrational effects at T 7^ K, however, lowers the Cu solu- 



data 



bility limit. It is not clear to us whether the reduction in solubility limits of Cu in bcc-Fe 

10 



calculated using the PM potential including the vibrational effects is real or an artefact of 
the employed potential. 




Present Calc. (with vib.) 
Present Calc. (no vib.) 
Pasianot 2007 (CO5.20_ATAT) 
Perez 2005 
Salje 1977 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Cu (at%) 



FIG. 4. Solubility limit of Cu in bcc-Fe. Present results are compared with experiments 



calculations performed with the ATAT package 



491] and 



501 1 using the PM potential denoted as CO5.20 



Nevertheless, after confirming the similar trend of total formation free energy of a single 
Cu-atom in bcc-Fe calculated with aforemetioned two different approaches namely, the har- 
monic approximation and the full thermodynamic integration, we further extended our stud- 
ies for pure Cu-clusters of larger sizes. The total formation free energy of pure Cu-clusters 
{Cum) increases with temperature within a range of 2 eV over the considered temperature 
range as shown in Fig. 5(b). This increase in Gj°*"^(Cum, T) means that formation or 
substitution of Cu-clusters in bcc-Fe matrix becomes somewhat more difficult at elevated 
temperature. 

Mixed vacancy-Cu clusters (viCum) as shown in Fig. 5(c) exhibit the effect of vibrational 
contribution (i.e. G^j^{viCumiV^T)) by reducing the values of their total formation free 
energy with temperature. The similar effect is also observed in the case of pure vacancy 
clusters (f/). Moreover, the number of vacancies present in viCum clusters seems to affect the 
formation energy of ViCu^ clusters (i.e. EfiyiCum)) prominently, and further influences the 
clustering of their formation energy, EfiviCum) (cf. Fig. 5 in 17j) and total formation free 



energy, Gj- {viCum, V, T). The vibrational contribution to the total formation free energy 



oi viCum clusters, Gf{viCum,V,T) appears to affect the trend of Gj- {viCum,V,T) as a 



function of temperature. The reduced values of G^°^°'''{viGum, V, T) at higher temperatures 
suggest the increased tendency of formation of mixed viGum clusters and such mixed clusters 



11 
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FIG. 5. Total formation free energy (eV) of (a) pure vacancy (vi) clusters (b) pure Cu (Cum) 
clusters and (c) mixed vacancy-Cu (viCum) clusters in bcc-Fe. 

are easier to form compared to pure vi clusters (cf. Fig. 5(a), Fig. 5(c)). In fact, in all cases 
discussed above, EjiyiCum) primarily determines the values of G^''^°'\viCum,V,T). The 
contribution from GJ^{viCum,V,T) further modifies these values for temperature effects 
and also determines the trend of G^°^"'''{viCum, V, T) with temperatures. 

C. Binding free energy of the clusters 

Figure 6(a) shows that the absolute value of total free binding energy for pure Vi clusters, 
G'™^'(f/, y, T) decreases with temperature. This decrease is associated to the increased 
vibrational contribution. The decreased absolute values of G]j^f{vi,V,T) indicate that the 
vibrational contribution weakens the bonding of pure vi clusters at elevated temperatures. In 
other words, pure vacancy clusters are less stable at elevated temperatures. At an operating 

12 



temperature of a nuclear reactor close to 600 K, the absolute values of total binding free 
energy are about 6% smaller than their corresponding values at K. 

Pure Cum clusters exhibit an increase in absolute values of both G'"f^^^{viCum, V, T) (not 
shown here) and G]j^f{viCum,V,T) as shown in Fig. 6(b) with increasing temperature 
which indicates that vibrational contribution strengthens the bonding of pure Cum clusters 
at elevated temperatures. This increase also means that in constrast with pure vacancy 
clusters, pure Cu clusters are more stable at higher temperatures. The absolute values of 
total binding free energy at about 600 K, are observed to increase by about 3-6% compared 
to EundiyiCum) at K. In the present study the absolute value of total free binding energy 
of a small Cuiq cluster, G^^fi^CuiQ, V, T) calculated at 800 K shows a good agreement with 



the first principle calculations reported at 773 K [20|. 

As shown in Fig. 6(c) the absolute values of G'[°^f{viCum,V,T) of all mixed vacancy- 
copper clusters considered here decrease with increasing temperature. At 1200 K the de- 
crease ranges from ~ 1 eV for viqCuiq clusters to ~ 12 eV for V8oCu2o- The V2QCum, v^oGum 
and VQoCum clusters exhibit a decrease of ~ 3-5 eV, ~ 6-8 eV and ~ 8-10 eV, respectively 
at 1200 K, where m = 20, 40, 60. This reduction is primarily attributed to the increase in 
Gll'^^{viCum, V, T) (not shown here) and appears to follow the trend of the total binding free 
energy of pure vi clusters with increase in temperature. In other words, the decrease in the 
absolute values of total binding free energy as a function of temperature indicates the higher 
probablity of dissociation of these mixed clusters into their constituents i.e. / monovacancies 
and m single Cu-atoms. At about 600 K, the absolute values of total binding free energy 
of viCum clusters is observed to reduce by minimum of ~ 2.5% for viqCuiq to maximum of 
about 6% for vsoGu2o in comparison with their corresponding values at K. The V2oGum, 
v^QCum and v^oCum clusters show a reduction of ~ 5-3%, ~ 6-4% and ~ 6-5%, respectively 
at 600 K in the absolute values of their Gj°l^f{viGum, V, T) compared to their corresponding 
values at K. 

In order to estimate the monomer binding energy for arbitrary pure cluster sizes, Al- 



Motasem et al. 



171] have fitted the dependence of total binding energy at K of pure 



clusters on cluster size to the following analytical relations: 



Eund{vi) ^ a f'^ + bl + c 

Eund{Gum) ^ dw?/^ + em + f (15) 
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FIG. 6. Total binding free energy (eV) of (a) pure vacancy (vi) clusters (b) pure Cu (Cum) clusters 
and (c) mixed vacancy-Cu (viCum) clusters in bcc-Fe. 

These fits, based on the classical capillary model, are valid for l,m > 2. The values of the 
fit parameters a (= 2.80595 eV), c (= -1.53677 eV), d (= 0.59667 eV) and / (= -0.60187 
eV) are determined by setting b and e to the negative formation free energies at K of the 
corresponding monomers, i.e. b = - Ef{vi) = -1.71 eV and e = - Ef{Cui) = -0.4369 eV. The 
main idea to determine these analytical fits is to provide sufficiently correct and compact 
description of the total and the monomer binding energy for further use in rate theory and 
object kinetic Monte Carlo simulations. The analytical relations in Eqs. ( |T5|) are chosen in 
such a way that the analytical form of their derivatives, i.e. 

E,{Cu, Cu^) ^ ^-^^-^^(^^^-) = Id m-V3 + e (16) 
corresponds to the relation for the monomer binding free energy used in conventinal rate 
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theory (cf. Eq.(4) in ref. jsL for OK). The coefficients a and d are related to the quantity 
^V^/^ {36n)^^^ given in ref. 5|, where 7 denotes the cluster- matrix interface energy and Vat 
the atomic volume in the bcc-Fe. 

Following the concept of Al-Motasem et al., we include here the temperature effects 
taking vibrational part of total binding free energy of pure clusters into consideration. The 
analytical expressions given by Eq. f llSp are then modified for temperature effects. They are 
re-written as 

^^(z;;, V, T) ^ a'(T) + b'{T) I + c'{T) 

Gl°"'\Cum, V, T) ^ d'{T) m2/3 + e'(T) m + /'(T) (17) 

The modified fit parameters are as follows: 

a'{T) = 2.76602 - 2.69198 x 10"^ T 

h'{T) = -G^^'^^ivi, V, T) = -1.70968 + 1.70404 x 10"^ T 

c\T) = -1.18921 + 1.07434 x lO"'^ T - 3.56495 x 10"^ 

d'{T) = 0.56749 - 4.38711 x 10"^ T 

e'(T) = -G^^^^^Cui, V,T) = -0.4382 - 4.0428 x 10"^ T 

f'{T) = -0.47811 + 1.82043 x 10~^ T - 4.02524 x 10^* T^. 

The monomer binding free energy at T 7^ and zero pressure defined in Eq. (IT^ can then 
be determined by treating Eq. f lT7|l with the similar analytical approximations introduced 
in Eq. ([16]). 



IV. SUMMARY AND CONCLUSION 

In this paper, we report our work on vibrational contributions to the thermodynamics 
of nano-sized vacancy-copper clusters in bcc-iron. The most recent many-body empirical 
potential developed by Pasianot and Malerba is employed. Vibrational effects are calculated 
under the harmonic approximation. We summarize our results as follows: 

1. The vibrational free energy of bcc-Fe and fcc-Cu compares well with SGTE database 
and data from the other empirical potentials. 
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2. The total formation free energy of a monovacancy in bcc-Fe shows good agreement with 



the work of Mendelev et ah 3J]. The reducing values of the total formation free energy 
for a monovacancy and pure vacancy clusters associated with increased vibrational 
effects at high temperatures exhibit an enhanced activity of vancancy formation in 
bcc-Fe at elevated temperatures. 

3. The total formation free energy of a single Cu atom in bcc-Fe contradicts the predic- 
tions of the first principle calculations and it is not clear to us whether the observed 
results is real or an artefact of the PM potential. The increase in the total formation 
free energy with temperature associated with increased vibrational contribution for a 
single Cu atom and pure Cu clusters indicates the increased difficulty of substitution 
of Cu atoms in bcc-Fe. 

4. The total formation free energy of mixed vacancy-copper clusters decreases with tem- 
perature due to increased vibrational contributions following a similar trend of pure 
vacancy clusters. For the number of vacancies and copper atoms considered in mixed 
clusters in the present work, vacancies appear to infiuence significantly the grouping 
of the total formation free energy of the mixed clusters. 

5. The observed decrease in the absolute values of total binding free energy of pure 
vacancy clusters with increasing temperature indicates a higher probability of the 
dissociation of pure vacancy clusters into individual vancancies. The absolute values 
of the total binding free energy of pure vacancy clusters show the decrease of about 
6% at an operating temperature of a nuclear reactor of about 600 K as compared to 
their corresponding values at K. 

6. The observed increase in the absolute values of the total binding free energy of pure 
Cu clusters on the other hand suggests a lower probability of the dissociation of the 
pure Cu clusters into individual components at increased temperatures. The absolute 
values of the total binding free energy of pure Cu clusters show the increase of about 
3-6% at about 600 K as compared to their corresponding values at K. 

7. The absolute values of the total binding free energy of mixed defect clusters follow the 
trend of pure vacancy clusters and show a decrease with increasing temperature. In 
case of mixed clusters with considered number of vacancies and copper atoms here, 
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vacancies overbalance vibrational contributions from copper atoms and exhibit an 
overall decrease in the absolute values of total binding free energy suggesting the higher 
probability of dissociation of mixed defects clusters into their individual constituents 
at higher temperatures. The absolute values of the total binding free energy of mixed 
clusters show a decrease from a minimum of ~ 2.5% to a maximum of ~ 6% depending 
on the number of vacancies present in the clusters. 

8. Finally, the temperature dependence is included in the analytical fits in order to pro- 
vide sufficiently correct and compact description of the monomer binding energy for 
further use in the rate theory and the object kinetic Monte Carlo simulations. 
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